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ABSTRACT 


We investigate the Sun-Earth dynamics of a set of eight well observed solar 
coronal mass ejections (CMEs) using data from the STEREO spacecraft. We 
seek to quantify the extent to which momentum coupling between these CMEs 
and the ambient solar wind (i.e., the aerodynamic drag) influences their dynam¬ 
ics. To this end, we use results from a 3D flux rope model fit to the CME data. 
We find that solar wind aerodynamic drag adequately accounts for the dynamics 
of the fastest CME in our sample. For the relatively slower CMEs, we find that 
drag-based models initiated below heliocentric distances ranging from 15 to 50 
Rq cannot account for the observed CME trajectories. This is at variance with 
the general perception that the dynamics of slow CMEs are influenced primarily 
by solar wind drag from a few Rq onwards. Several slow CMEs propagate at 
roughly constant speeds above 15-50 R e . Drag-based models initiated above 
these heights therefore require negligible aerodynamic drag to explain their ob¬ 
served trajectories. 


Subject headings: Sun: Corona — Sun: Coronal mass ejections (CMEs)—Solar wind 


1. Introduction 


Coronal Mass Ejections (CMEs) are hot blobs of plasma and magnetic field that 
occasionally erupt from the solar corona. Some CMEs are directed towards the Earth 
and are primary drivers of the majority of so-called “space weather” disturbances that 
have profound consequences on a variety of technologies that we use in everyday life. 
Space weather predictions are therefore crucially dependent on accurate estimates of the 
Sun-Earth travel time of Earth-directed CMEs and their speed upon arrival at the Earth. 
Although we have some understanding of the factors influencing CME propagation through 
the inner heliosphere, several crucial aspects such as the origin of the driving and drag 
forces are still somewhat unclear. Needless to say, a thorough understanding of the physical 
bases of CME propagation is crucial to enable reliable space weather predictions. 

Research on CMEs is often divided into “initiation” and “propagation” issues. 
Initiation issues deal with the initial eruption stages of a CME, while propagation issues 
concern forces at work during its subsequent propagation through the solar corona and 
the heliosphere. The division is somewhat artificial, but one can think of propagation 
issues as the ones that deal with the post-initiation, “residual” acceleration stage of CMEs 
(e.g., Zhang & Dere 2006; Subramanian & Vourlidas 2007; Gopalswamy 2013). Based on 
the temporal evolution of the kinetic and potential energies of a sample of well observed, 
plane-of-sky CMEs, Subramanian & Vourlidas (2007) showed that they are acted upon 
by a driving force throughout the Large Angle and Spectrometric Coronagaph (LASCO; 
Brueckner et al. 1995) field of view (FOV). Around 17 % of CMEs in the extensive 
LASCO dataset are observed to accelerate out to 30 R Q (St.Cyr et al. 2000). We will 
concern ourselves here with the propagation stage of CMEs. During the propagation stage, 
CMEs are thought to be subject to Lorentz self-forces close to the Sun ( which may be a 
continuation of the initiation phase) and to aerodynamic drag due to momentum coupling 
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with the ambient solar wind farther out. The interplay between these forces determines the 
dynamics of a CME as it propagates outwards from the Sun. It is generally accepted that: 

1. very fast CMEs (generally defined as ones with speeds in excess of 1000 kms -1 near 
the Sun) experience all their acceleration (due to Lorentz self-forces) very close to 
the Sun, and are subsequently subjected mostly to deceleration due to aerodynamic 
drag from the ambient solar wind; i.e., they are dragged down by the solar wind via 
aerodynamic drag. 

2. the dynamics of slow CMEs (generally defined as CMEs with speeds comparable to 
the solar wind speed near the Sun, of the order of a few hundred kms -1 ) are governed 
primarily by solar wind aerodynamic drag. While Lorentz self-forces are expected to 
play some role, it is often thought that slow CMEs are primarily dragged up by the 
solar wind. 

Our understanding of how Lorentz self-forces operate within CMEs is not very robust, 
although we know that they can accelerate, as well as decelerate a CME. While Lorentz 
self-forces must involve misaligned currents and magnetic fields, it has only recently been 
demonstrated that flux rope CMEs need to be substantially non force-free in order to 
account for Lorentz self-force driving (Subramanian et al. 2014). There is observational 
evidence to argue that the magnetic energy contained in flux rope CMEs is adequate to 
account for the mechanical energy required to propagate them from the Sun to the Earth 
(e.g., Manoharan et al. 2004; Manoharan 2006; Subramanian & Vourlidas 2007). However, 
the actual nature of Lorentz self-forces has not yet been well understood; some authors 
incorporate it in the eruption process (Isenberg & Forbes 2007; Kliem et al. 2014), while 
some treatments (e.g., Chen 1996) appeal to poloidal flux injection as a proxy for Lorentz 
self-forces, and use soft X-ray flare profiles as guides to tailor the time profile of poloidal 
flux injection (Chen & Kunkel 2010). 
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On the other hand, the idea of momentum coupling with the ambient solar wind 
was motivated by the findings of Gopalswamy et al. (2000), Manoharan (2006), Maloney, 
Gallagher & McAteer (2009) and others who found that CMEs launched with speeds 
exceeding that of the solar wind were decelerated in the interplanetary medium, while those 
launched with speeds slower than the solar wind were sped up, so that all CMEs tend to 
equilibrate to the speed of the solar wind. Using a novel method of identifying an “average” 
CME from LASCO data over a one year period during the approach to solar maximum, 
Lewis & Simnett (2002) arrived at similar conclusions based on the similarities between the 
velocity profile of the CME and the ambient solar wind. Such observations have spurred 
considerable research on drag-based models (e.g., Cargill 2004; Vrsnak 2006; Howard et 
al. 2007; Borgazzi et al. 2009; Byrne et al. 2010; Maloney & Gallagher 2010;Vrsnak et al. 
2010; Temmer et al. 2011; Lugaz & Kitner 2013; Mishra & Srivastava 2013; Dolei et al. 
2014; Vrsnak et al. 2013 ; Iju,Tokumaru & Fujiki 2014; Temmer & Nitta 2015). Most of 
these investigations employ an empirical drag coefficient Cx> in order to characterize the 
interaction between the CME and the solar wind. An exception is Cargill (2004) where Cb 
is motivated from the results of 2.5D MHD simulations. Subramanian, Lara & Borgazzi 
(2012) (SLB2012 hereafter) have obtained a physical prescription for Cb from a calculation 
of collisionless viscosity in the solar wind. The Sun-Earth travel time and near-Earth 
CME speed predicted by a model using this prescription is found to agree very well with 
observations of fast CMEs (SLB2012). 

Using a large sample of CMEs observed in LASCO aboard the SOlar and Heliospheric 
Observatory mission (SOHO; Domingo, Fleck & Poland 1995) mission, Michalek, 
Gopalswamy & Yashiro (2015) have recently characterized three phases of propagation. 
The first phase is the initial acceleration phase (presumably) due to Lorentz self-forces, 
followed by the second phase where a balance between the Lorentz and drag forces acting 
on the CMEs is achieved. The third phase is when the CME primarily decelerates due to 


aerodynamic drag at the outer edge of LASCO FOV. Our approach is similar, but we only 
seek to differentiate when CMEs are drag-dominated and when they are not. We note that 
the solar wind drag can serve to decelerate, as well as accelerate CMEs. 


2. Data and Models 
2.1. CME Event Sample 

We use a set of eight extensively studied, Earth-impacting CMEs during the rising 
phase of solar cycle 24 between March 2010 and June 2011 (Colaninno 2012). These CMEs 
are observed by the SOHO LASCO, Sun-Earth Connection Coronal and Heliospheric 
Investigation (SECCHI; Howard et al. 2008) coronagraphs and the Heliospheric Imagers 
on board the Solar TErrestrial RElations Observatory mission (STEREO; Kaiser et al. 
2008). Near-Earth parameters for these CMEs are obtained via in-situ measurements 
from the WIND spacecraft (http://omniweb.gsfc.nasa.gov). The data set is described 
comprehensively in Colaninno (2012) and Colaninno, Vourlidas & Wu (2013), who have 
used the Graduated Cylindrical Shell Model (GCS) (Thernisien, Howard & Vourlidas 2006; 
Thernisien, Vourlidas & Howard 2009) to obtain a 3D reconstruction of the Sun-Earth 
trajectory for each of these CMEs. This represents one of the best samples of CMEs for 
which detailed 3D information is available from the Sun to near the Earth. The ninth 
CME in the dataset of Colaninno (2012) involves CME-CME interactions, which cannot be 
accommodated by the simple, drag-based model we employ. We therefore omit that event 
and study only eight CMEs from the sample. 

Table 1 contains salient details for each of these CMEs. The CMEs are hereafter referred 
to using the serial number assigned in this table. For each CME, Colaninno (2012) gives 
the detailed 3D height-time trajectory and flux rope morphology from the Sun to near the 
Earth. In addition, we use some other parameters which we describe herein. The quantity 
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ho denotes the height at the first measurement and hf denotes the height at the final 
measurement. The quantity Vq is the initial speed of the CME. We note that CMEs 1, 4, 6 
and 8 start out fairly slow, with initial speeds < 220 km s _1 , while CME 3, 5 and 7 start 
out relatively faster (>400 km s” 1 ). CME 7 was preceded by another CME, thus the solar 
wind characteristics are expected to be different from the quiescent state. CME 2 starts out 
with an initial speed of 916 km s -1 and is the fastest event in this sample. 


2.2. Drag Model 


Momentum coupling between the CME and ambient solar wind is typically addressed 
using the following equation; 


h'drag ^Hcme 


d.V„ 


dt 


2 k V) A erne ( V C jnp Ik, 


keme V sw 


( 1 ) 


where Fd rag is the drag force due to momentum coupling between the CME and solar wind, 
m cme is CME mass, V cme is CME speed, Cd is the drag coefficient, A cme is CME area, n is 
the solar wind density, m p denotes the proton mass and V sw is the ambient solar wind speed. 
The dimensionless drag parameter (Cd) expresses the strength of the momentum coupling 
between the CME and the solar wind. In fluid dynamics texts (e.g., Landau & Lifshitz 
1987), an equation of motion like Equation (1) is used to describe high Reynolds number 
flows past a solid body, where the bulk of the flow (which can be regarded as a largely 
potential flow) is separated from the body by a turbulent boundary layer. The solid body 
boundary conditions require that the entire velocity (both the tangential, as well as the 
normal component) vanish at the boundary. The turbulence in the boundary layer (which 
is rather thin, owing to the large Reynolds number) is set up because the velocity needs to 
transition from its bulk value to zero on the boundary over a very short distance; nonlinear 
terms in the Navier-Stokes equation assume importance and give rise to turbulence. 
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We note that Equation 1 is also often expressed in the following form (e.g., Cargill 
2004; Vrsnak et ah 2010): 



( 2 ) 


where, 7 (cm x ) is given by 


nm n A , 


( 3 ) 


The cross-sectional area of CME is calculated using A cme (R cme ) = tt R cme W cme , where 
Rome is the minor radius and W cme is the major radius of the elliptical cross section in the 
x-z plane of a flux-rope CME (Thernisien 2011). For each CME in our sample, the GCS 


fitting yields the major radius, the ratio of the minor to major radius and the half angle 


which can be used to calculate R cme and W cme using equations 27, 28 and 29 of Thernisien 
( 2011 ). 


In this paper we will use the SLB2012 prescription for Cd, as well as prescriptions 
where the value of Cd is maintained at a constant value. For the sake of completeness, 
we briefly recapitulate the essentials of the SLB2012 viscosity and Cd prescription. The 
viscosity (u sw ) in the ambient collisionless solar wind is due to resonant scattering of solar 
wind protons with an Alfven wave spectrum, and can be expressed in analogy with a fluid 
expression for turbulent viscosity (Verma 1996): 



( 4 ) 


where u rms is the rms speed of solar wind protons and A is the mean free path, which we 
take to be the ion inertial length, (Leamon et al. 1999, 2000; Smith et al. 2001; Bruno & 
Trenchi 2014) and is given by 


c 


A 


r\j 


228 n x / 2 km . 


( 5 ) 




9 


In Equation (5), v a is the Alfven speed, fh is the ion cyclotron frequency, uj p is the ion 
plasma frequency and n is the ambient solar wind density. We note that this prescription 
is similar to that used by Coles & Harmon (1989); but their prescription is a factor of 3 
larger. The treatement of SLB2012 uses the Coles & Harmon (1989) prescription for A, 
whereas we take it to be the ion inertial length. The CME Reynolds number is defined as 

Re = {VcME ~ Vm) Rcme . (6) 

VSW 

where the quantity v sw (cm 2 s" 1 ) is the solar wind viscosity calculated using the prescription 
given in SLB2012. The drag parameter Cd is related to Reynolds number using a fit to the 
data from Achenbach (1972) for the drag on a sphere at high Reynolds numbers: 

C D = 0.148 - 4.3 x 10 4 Re" 1 + 9.8 x 10" 9 Re . (7) 

The Cd prescription given by Equation (7) was experimentally determined for subsonic, 
high Reynolds number flow past a solid metal sphere (Achenbach 1972). It is worth 
commenting on its applicability to our situation. Firstly, as noted earlier, the basic drag 
force law (Equation 1) that is used extensively in this held is really a high Reynolds number, 
solid body law - two of the assumptions inherent in Equation (7) are thus consistent with 
it. Furthermore, it has been shown that the total (magnetic + particle) pressure exhibits 
a substantial jump across a typical magnetic cloud boundary, suggesting that they are 
over-pressured structures (Russell, Shinde & Jian 2005; Jian et ah 2006) that will not 
deform in response to tangential stresses, much like a solid body. Achenbach’s drag formula 
has also been verified by modern detached-eddy simulations of high Reynolds number hows 
in the super-critical regime (e.g., Const ant inescu & Squires 2004), which is the regime we 
are interested in. Another aspect to address is the fact that Equation (7) is derived from 
experimental results for subsonic how, whereas the motion of CMEs through the solar wind 
is certainly supersonic. According to the Morokovin hypothesis, which has been verihed 
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extensively via numerical simulations (e.g., Duan, Beekman & Martin 2011) a subsonic 
turbulent drag law is valid even for supersonic flows as long as the fluctuations in the 
turbulent boundary layer are incompressible, or subsonic. The turbulent fluctuations in the 
ambient solar wind are certainly incompressible to a very good approximation (e.g., Shaikh 
& Zank 2010). This is also evident from the small values of the density fluctuation A n/n 
in the ambient solar wind (e.g., Bisoi et al. 2014). The values for the fluctuations in the 
magnetic field A B/B in the turbulent sheath between the shock and the CME are also 
as small as 10% (Arunbabu et al. 2013). Since A B/B is generally representative of A n/n 
(e.g., Spangler 2002), it follows that the turbulent density fluctuations in the sheath region 
are also fairly incompressible. In summary, the SLB2012 Cb prescription, as described by 
Equations (4), (5), (6) and (7) is quite appropriate to our situation. 

However, as mentioned earlier, we use the SLB2012 Cb prescription as well as constant 
values for Cd in the results described below. 


2.3. Models for solar wind density, speed and CME mass 


We next describe the data-based models we use for computing the CME mass (m cme ), 
solar wind speed {V sw ) and solar wind density (n) in Equation 1. The quantity n W i n d in 
Table 1 denotes the proton number density at the Earth, approximately one to two days in 
advance of the CME and the shock arrival. Assuming that the proton and electron number 
densities are equal, we use a version of the model given by Leblanc et al. (1998) in order 
to extrapolate the proton number density Sunwards from the Earth. The electron number 
density at a given heliocentric distance is given by: 


n(R ) = 


( W"wind \ 

~tYJ [ 


3.3 x 10 5 ir 2 + 4.1 x 10 °ir 4 + 8 x 10 'R 


6 D—4 


W D—6 


cm 


-3 


( 8 ) 
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where R is the heliocentric distance in solar radii. The original model of Leblanc et 
ah (1998) assumes that the number density at 1 AU is 7.2 cm -3 . The factor n vnnr ]j7.2 
ensures that the number density given by Equation 8 matches the proton number density 
(■ n W ind ) measured in-situ by the WIND spacecraft near the Earth. The quantity n(R) 
thus represents the number density of the ambient solar wind into which a given CME 
propagates. 

It is well known that the ambient solar wind speed is an important factor influencing 
CME propagation (eg., Tennner et al. 2011). The quantity v W i n d in Table 1 denotes the 
near-Earth solar wind speed observed in-situ by the WIND spacecraft, in advance of the 
CME arrival. In order to obtain the solar wind speed V SW (R) as a function of heliocentric 
distance (R) along the Sun-Earth line, we use the solar wind model of Sheeley et al. (1997, 
1999): 

V,l(R) = vl ind [1 - e -^] . (9) 

The quantity ro, which denotes the heliocentric distance where the solar wind speed is zero, 
is taken to be 1.5 R 0 . The quantity r a , which is the e-folding distance over which the 
asymptotic speed v win d is reached, is taken to be 50 R e . We have veriEed that our results 
are not very sensitive to the precise value of r a ; this will be described later in this paper. 

To calculate the CME mass, we used the Thomson scattering geometry from Billings 
(1966), which relates the observed brightness to coronal electron density. We calculated the 
mass in the COR2 FOV assuming that all the mass of the CME is concentrated in a plane 
at the longitude derived from the GCS model. This assumption is a first order estimate of 
the 3D mass consistent with the method of Colaninno & Vourlidas (2009). To correct for 
the effects of the occulter on the observed 3D mass, we used the prescription of Bein et al. 
(2013) who suggest the following functional form in the coronagraphic FOV: 


m cme (R) = m 0 + Am (R - h 

occ) • 


( 10 ) 
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In this equation, mo is the “true” mass of the CME when it is first visible above the 
occulter. The quantity Am is the real mass increase with height and h occ is the size (radius) 
of the occulter in projection (Table 1). These terms in Equation 10 are derived from fitting 
the observed mass increase as the CME emerges from behind the occulter. Thus we can 
estimate the un-occulted 3D mass of the CME in the COR2 FOV with Equation 10. This 
mass formulation is valid only within the COR2 FOV, therefore, we assume the mass to be 
constant beyond 15 R 0 . 

We have also added the virtual mass (e.g. Cargill et al. 1996; Cargill 2004) to m crne in our 
calculations. We note that the addition of the virtual mass does not change the results in 
any significant manner. 


3. Model predictions vs data 

To recapitulate, Equation 1 is a simplified, one-dimensional equation of motion 
describing the dynamics of a CME that is governed solely by aerodynamic drag due to 
momentum coupling with the ambient solar wind. Other effects such as Lorentz self-force 
driving are not included in this treatment. Due to momentum coupling between the CME 
and the solar wind, CMEs that travel slower than the solar wind would thus be accelerated, 
while those travelling faster than the solar wind would be decelerated. As described above, 
we have used realistic observational values for all parameters. The CME speed supplied 
to the model as an initial condition is estimated from the observed height-time data 
using either a multipolynomial fit as in Colaninno (2012) or via numerical differentiation. 
We compare the prediction of the CME trajectory from Equation 1 with the observed 
height-time trajectory obtained via 3D GCS flux rope fitting for each of the CMEs in 
Table 1. Equation 1 yields the speed of a CME (V cme ) as a function of time from the start 
of the event. The output of Equation 1 is integrated to obtain a height-time trajectory 
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predicted by the model. We compare the predicted height-time trajectory with the observed 
height-time data for each CME in the next subsection. Other relevant parameters dealing 
with the comparison between the observations and corresponding models are discussed in 
$ 3.1. 


3.1. Comparision of height-time profiles 

The observed height-time data together with the model predictions are shown in 
Figures 1 and 2. The results for CMEs 1, 2, 3 and 4 of Table 1 in Figure 1 and those for 
CMEs 5, 6, 7 and 8 are shown in Figure 2. Observed height-time data points are denoted 
by diamonds, while the lines show the model predictions. The dash-dotted line (red) shows 
the prediction of the model (Equation 1) when it is initiated with quantities corresponding 
to the first data point. The solid (blue) line, on the other hand, shows the model prediction 
when it is initiated with these quantities corresponding to a later time; in other words, the 
model CME is initiated from a height larger than the first observation. Other details about 
the model predictions for each CME are discussed in the following sections. 


3.1.1. CME 2 

We note that the dash-dotted line which is computed using the SLB2012 Cd 
prescription agrees reasonably with the data points for the fastest CME in our sample, 
CME 2, which has a starting speed of 916 kms -1 . This is expected, since fast CMEs, which 
are primarily decelerating, are dominated by solar wind aerodynamic drag, which is the 
only force that is included in Equation 1. 

Although the height-time plot (Figure 1) looks roughly like a straight line, and might 
suggest that the CME speed is nearly constant, we note that this is not so; CME 2 
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decelerates from ~ 916 kms 1 to ~ 715 kms 1 . 

The modeled CME is somewhat faster than the observed one; the observed Sun-Earth 
travel time for CME 2 was ~ 60 hours and the model predicts that the CME arrives at the 
Earth around 8 hours earlier. This represents an error of 14%. Ah which is the difference 
between the last measured height from observations and the final modeled height, was -31.1 
R 0 for this CME. 

Our time of arrival (ToA) errors are comparable to the values of « 10-12 hours obtained 
from MHD models (Mays et al. 2015) and drag-based models (Shi et al. 2015). Although we 
have used observational data as much as possible for the model fitting, there is some room 
for uncertainty in the quantity r a . We find that a decrease (increase) in r a of 50% results in 
a 3.7% decrease (increase) in the predicted CME travel time, ft is also possible that errors 
in the GCS flux rope fitting procedure can lead to errors in the measured cross-sectional 
area of the CME. We find that a 50 % decrease in the CME area results in a 6 % decrease 
in the CME travel time. 

The range of Cd values predicted by STB2012 Cd prescription for CME 2 is 0.6-1.4. Model 
predictions using a constant 0.6 < Cd < 2.0 are reasonably similar to the dash-dotted 
(red) line of CME 2 in Figure 1 which uses the STB2012 Cd- On the other hand, constant 
Cd models with values for Cd that lie substantially outside of 0.6 < Cd < 2.0 perform 
poorly. This is evident from the examples shown in Figure 1; model solutions with constant 
Cds of 0.1 (green dash-dotted line) and 5 (brown dash dotted line) disagree considerably 
with the observed height-time data for CME 2. For the Cd = 0.1 model, the difference 
between the predicted and the observed ToA at 180 R 0 is ~ -8 hours. For Cd = 5.0, it is 




14 hours, and only -1.8 hours for the SLB2012 Cd prescription. 
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3.1.2. Slower CMEs 

Next, we turn our attention to the somewhat slower CMEs in our sample (Table 1): 
CMEs 1, 3, 4, 5, 6, 7 and 8. For such CMEs, it is generally believed that Lorentz self-forces 
are dominant for the initial part of its trajectory, while aerodynamic drag takes over at 
larger heights (e.g., Michalek, Gopalswamy & Yashiro 2015). Some authors (e.g., Lewis 
& Simnett 2002) claim that representative slow CMEs are exclusively dragged up by the 
solar wind. Some authors (e.g., Mishra & Srivastava 2013) apply the drag-based model to 
explain observations of slow CMEs over its entire trajectory while some authors (e.g., Byrne 
et al. 2010; Carley, McAteer & Gallagher 2012; Gopalswamy 2013) think that solar wind 
aerodynamic drag becomes dominant for slow CMEs beyond a few solar radii. It is evident 
from Figures 1 & 2 that the dash-dotted line disagrees considerably with the data (denoted 
by diamonds) for all events except CME 2. The maximum height (at the last timestamp) 
predicted by the model (when initiated from the start) falls short of the final observed 
height (hf) for each CME (except CME 7, in which case it is much larger than hf). We find 
that the difference between the final height predicted by the model and the observed final 
height is large, ranging between 17-175 R 0 for CMEs 1, 3, 4, 5, 6, 7 & 8. This suggests 
that there is a large discrepancy in the predicted and observed trajectory when the model 
is initiated from low heights. In other words, momentum coupling with the ambient solar 
wind alone does not satisfactorily explain the dynamics of the CMEs other than CME 2. 

While the dash-dotted (red) lines use the SLB2012 Cd prescription, we also find that 
constant Cd models with 0.1 < Cx> < 5 cannot account for the data when initiated from 
the first observation. This points to the conclusion that CMEs 1, 3, 4, 5, 6, 7 and 8 cannot 
be considered to be drag-dominated (either drag-accelerated or decelerated) from the start. 
This conclusion is independent of the specific Cd model used (the SLB2012 Cd or a constant 
Cd). In order to find the height where the solar wind drag “takes over” for the relatively 
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slower CMEs we adopt the following strategy - we initiate the model at a larger height. 

The initial conditions supplied to the differential equation (Equation 1), such as the CME 
initial speed, mass, cross-sectional area, solar wind speed, background density are those 
appropriate to this height, which is larger than that at the first observed data point. We 
repeat this exercise with increasing initiation heights until there is a reasonable agreement 
between the modeled (solid line) and observed height-time data. We carry out this exercise 
using the SLB2012 Cd model as well as constant Cd models. As with CME 2, we find that 
the values for Cd predicted by the SLB2012 prescription are a good guide for the constant 
Cb model. In other words, only those constant Cd models with Cd values that are close 
to those predicted by SLB2012 prescription (last column in Table 2) agree reasonably well 
with the data. The initiation height beyond which the solid line agrees reasonably with the 
data is denoted by h 0 in Table 2, and the CME speed at h 0 is denoted by v 0 . At heights 
below h 0 , no drag model yields good agreement with observations. We have verified this 
by initiating the model at different heights ranging from the first observed height (ho) to 
ho, using the SLB2012 Cd prescription, as well as the constant Cd prescription. Since we 
do not include any forces other than solar wind drag, the only definitive statement we can 
make about CMEs 1, 3, 4, 5, 6, 7 and 8 is that they are not drag dominated when the 
model is initiated at any height below h 0 . 

Table 2 also lists for each CME the quantity hf which is the final height predicted by 
the model when it is initiated from ho (corresponding to the solid lines in Figures 1 and 2 
and dash-dotted line for CME 2). Ah is the difference between h/ and h/. 

For CMEs 1, 3, 4, 5 and 6, we have reconstructed 3D data only up to ~ 160 R & . In 
order to avoid empirical extrapolation, we therefore compare the model velocity predictions 
with the observed quantities at the last timestamp for all CMEs. Table 3 lists the quantities 
Vidp, which is the velocity observed at the last timestamp (at h/) for each event and V pre d, 
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which is the predicted velocity at the same instant. AC is the difference between V pre d and 
Vidp■ A positive AC means that the predicted CME velocity at last data point is larger 
than the observed one, while a negative AC means the converse. The ratio of AV to the 
observed velocity Cdp (at heliocentric height hf) is also given in Table 3. A V/Vd p ranges 
between 0.3% to 31% which indicates that the model performs well. Depending on how the 
ambient solar wind speed compares with v 0 , the solar wind drag decelerates/accelerates the 
CMEs. 

As with the velocity comparisons, in order to avoid empirical extrapolation, we 
compare the Time of Arrival (ToA) for each event at two heights: 130 R 0 and 180 R 0 . 
Only CMEs 2, 3, 7 and 8 attain heights > 180 R 0 ; therefore ToA comparisons at 180 R 0 
are made for only these four events. For the remaining events we make ToA comparisons at 
130 R 0 . Table 3 describes the quantity AT which is the difference between the predicted 
and observed ToA at these heights, when the model is initiated at h 0 (except for CME2, 
where the model is initiated at h 0 ). A positive AT means that the arrival time of the model 
CME is delayed with respect to the observed one, while a negative AT denotes its converse. 
The quantity AT/T denotes the ratio of AT to the observed CME transit time up to each 
height. It ranges between 1% and 5% for 130 R 0 and for 180 R 0 it ranges between 0.1% 
and 4%. 

An important point to be noted is that that most of the slow CMEs show little 
evolution in their speeds beyond h 0 . While CME 5 decelerates moderately (~ 13 km s -1 ) 
between 28 R e and 150 R e , CME 1 hardly decelerates. As before, constant Cd models 
with values of Cd that are close to the SLB2012 Cd predictions yield similar results. Since 
the deceleration is so small, even constant Cd models with Cd ~ 0 yield good results. 
Initiating drag-based models above h 0 thus does not constrain the drag-based models to 
any appreciable extent, and only reinforces the fact that most slow CMEs don’t accelerate 


or decelerate much beyond ho. 


4 . Discussion and Conclusions 

We worked with a dataset of eight well observed CMEs that were tracked by the SOHO 
LASCO and SECCHI coronagraphs and the HI imagers aboard the STEREO spacecraft 
(Colaninno 2012). These events were fitted with the GCS model, which yields the 3D 
height-time profile of the CME as it propagates through the heliosphere (Table 1). Despite 
the relatively small number of events, our sample covers a wide range of speeds: some 
events start out fairly slow (around a few hundred kms -1 ), some have intermediate speeds 
and one is relatively fast (starting speed of 916 kms -1 ), although it is certainly not among 
the fastest CMEs observed. We focused here on the effect of solar wind aerodynamic 
drag on CMEs using a simple, widely used ID model (Equation 1) that allows one to 
obtain the predicted height-time trajectory of a CME that is subject only to this drag 
(Lorentz self-forces are not included). It may be noted that solar wind drag can serve to 
decelerate, as well as accelerate a CME. We used two prescriptions for the drag parameter 
Cx>: the SLB2012 prescription, and a prescription where the value of Cx> is maintained at a 
constant value. For a given CME, we use observationally determined values for the CME 
cross-sectional area, mass, initial speed, ambient solar wind density and speed to solve 
Equation 1 and obtain the predicted height-time trajectory, which we compare with the 
observed one. In general, such a model is expected to work best for very fast CMEs, for 
which (accelerating) Lorentz self-forces are expected to cease very soon after initiation (e.g., 
Chen & Kunkel 2010), and the only operating force from then onwards is the (decelerating) 
solar wind aerodynamic drag. In keeping with this expectation, we find that the drag-based 
model performs well in describing the dynamics of the fastest CME in our sample, CME 2. 

Its generally thought that slow CMEs are “picked up by the solar wind” (i.e., 
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accelerated by aerodynamic drag) from as early as a few R®. Results from drag-based 
models reported in the most of the existing literature refer to initiation from the first 
observed timestamp, which typically correspond to only a few R®. For the slow CMEs 
in our sample, we find that the predictions of drag-based models do not agree with the 
observations when initiated at any heliocentric height below ho (ho for slow CMEs ranges 
from 15-50 R®). When initiated at heights < ho, the errors in predicting the final height 
for slow CMEs range from 12 % to 80 %. Our results thus suggest that solar wind drag is 
not the dominant influence on CME trajectories until 15-50 R®. It might be important to 
include other forces such as Lorentz forces, at least as far as 15 to 50 R®. Recent work by 
Zic, Vrsnak & Tcnnner (2015) also assumes that solar wind drag is effective only above 15 
R©. 

The speeds of several slow CMEs in our sample stay roughly constant for heights > h 0 . 
Although drag-based models initiated above h 0 seem to generally perform well, solar wind 
drag in fact “doesn’t have much to do” for these CMEs (since the CME speed evolution at 
heights > ho is negligible) and the drag parameter Cd is not well constrained. 
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CME 1 CME 2 






Fig. 1.— Height-time observations and model predictions for CMEs 1, 2, 3 and 4. The 
height-time data points are denoted by diamond symbols along with error bars for COR2 & 
HI data. The red dash-dotted line denotes the predictions of the model when it is initiated 
from the first observed data point, while the blue solid line denotes the model predictions 
when it is initiated from h 0 . In the case of CME 2, the green dash-dotted line indicates the 
model solutions when initiated from the start with a constant Cb of 0.1 while the brown 
dash-dotted line represents the model predictions using a constant Cb of 5.0. 
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CME 5 




CME 6 




Fig. 2.— Height-time plot for data and derived heights from the drag-only model for CMEs 
5, 6 ,7 and 8. The height-time data points are denoted by diamond symbols along with 
error bars for COR2 & HI data. The dash-dotted lines denoted the predictions of the model 
when it is initiated from the first observed data point, while the solid line denotes the model 
predictions when it is initiated from h 0 . 






Table 1. 


Observed CME and solar wind parameters 


No. 

Date 

ho 

(R©) 

h f 

(Re) 

vo 

(kms -1 ) 

V'wind 

(cm -3 ) 

Vyjind 

(kms -1 ) 

log (mo) 

h-occ 

(R-o) 

log(Am) 

1 

2010 Mar 19 - 2010 Mar 23 

3.5 

155.5 

162 

3.60 

380 

14.7 

3.65 

14.3 

2 

2010 Apr 03 - 2010 Apr 05 

5.5 

210.0 

916 

7.15 

470 

15.1 

3.46 

14.6 

3 

2010 Apr 08 - 2010 Apr 11 

2.9 

195.1 

468 

3.60 

440 

15.4 

3.55 

14.7 

4 

2010 Jun 16 - 2010 Jun 20 

5.7 

163.3 

193 

3.50 

500 

14.9 

4.37 

14.1 

5 

2010 Sep 11 - 2010 Sep 14 

4.0 

150.3 

444 

4.00 

320 

15.2 

3.84 

14.7 

6 

2010 Oct 26 - 2010 Oct 31 

5.3 

131.5 

215 

3.80 

350 

15.2 

5.66 

14.7 

7 

2011 Feb 15 - 2011 Feb 18 

4.4 

188.6 

832 

2.50 

440 

15.6 

3.85 

13.9 

8 

2011 Mar 25 - 2011 Mar 29 

4.8 

211.0 

47 

3.00 

360 

15.4 

4.43 

14.3 


Note. — Column 1: Serial number; Column 2: Observed date range for the event; Column 3: Initial observed height; 
Column 4: Final measured height; Column 5: Initial CME speed; Column 6: Observed proton number density at 1 AU; 
Column 7: Observed solar wind speed at 1 AU; Column 8: log of the “true” CME mass at the first measurement; Column 
9: Occulation height; Column 10: log of mass increase per height 
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Tabic 2. CME model parameters and SLB2012 Cd 








No. 

h 0 

v 0 

hf 

Ah 

SLB2012 C D 


(Re) 

(kms -1 ) 

(Re) 

(Re) 


1 

21.9 

383 

148.3 

7.1 

0.1-0.3 

2 

5.5 

916 

241.1 

-31.1 

0.6-1.4 

3 

19.7 

506 

191.8 

3.3 

0.2-0.4 

4 

15.2 

437 

158.4 

4.8 

0.2-0.3 

5 

27.7 

490 

149.1 

1.2 

0.4-0.6 

6 

20.1 

445 

141.7 

-10.3 

0.25-0.5 

7 

39.7 

530 

190.7 

-2.1 

0.25-0.4 

8 

46.5 

456 

212.0 

-1.0 

0.25-0.35 


Note. — Column 1: CME number corresponding to 
Table 1; Column 2: Height at which model is initiated; 
Column 3: CME speed at height ho; Column 4: Final 
height at last observed timestamp when the model is ini¬ 
tiated from h 0 ; Column 5: Difference between hf and hf] 
Column 6: Range of Cn values predicted by SLB2012 Cn 
prescription. 
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Table 3. CME arrival time and speeds : observations vs predictions 


No. 

Vid P 

(km s -1 ) 

Vpred 

{km s -1 ) 

AV 

(km s _1 ) 

AV/V ldp 

% 

130 R 0 

AT A T/T 

hr % 

180 R 0 

AT A T/T 

hr % 

1 

439 

382 

-57 

13 

0.85 

1 



2 

639 

731 

92 

14 



-1.8 

4 

3 

501 

503 

2 

0.3 



0.07 

0.09 

4 

501 

483 

-18 

3 

2.5 

4.5 



5 

564 

477 

-87 

15 

-1.98 

3.8 



6 

457 

439 

-18 

4 

-0.75 

1.27 



7 

595 

526 

-69 

12 



0.12 

0.2 

8 

342 

451 

109 

31 



-0.39 

0.46 


Note. — Column 1: CME number correspoing to Table 1; Column 2: Observed 
CME velocity at last observed timestamp; Column 3: Predicted model velocity at the 
last timestamp when the model is initiated from h 0 (except CME 2 in which case model 
starts from ho); Column 4: AV is the difference between V pre d and V[d p ] Column 5: Ratio 
of AV and Vi dp ; Column 6: Difference between the Predicted and Observed ToA at 130 
R 0 ; Column 7: Ratio of AT and the observed ToA at 130 R 0 ; Column 8: Difference 
between the Predicted and Observed ToA at 180 R 0 ; Column 9: Ratio of AT and the 
observed ToA at 180 R 0 . 






